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Abstract. 

We investigate the influence of a Bose-Einstein condensate (BEC) on the 
properties of immersed impurity atoms, which are trapped in an optical lattice. 
Assuming a weak coupling of the impurity atoms to the BEC, we derive a 
quantum master equation for the lattice system. In the special case of fixed 
impurities with two internal states the atoms represent a quantum register and 
the quantum master equation reproduces the exact evolution of the qubits. We 
characterise the qubit dephasing which is caused by the interspecies coupling 
and show that the effect of sub- and superdecoherence is observable for realistic 
experimental parameters. Furthermore, the BEC phonons mediate an attractive 
interaction between the impurities, which has an important impact on their 
spatial distribution. If the lattice atoms are allowed to move, there occurs a 
sharp transition with the impurities aggregating in a macroscopic cluster at 
experimentally achievable temperatures. We also investigate the impact of the 
BEC on the transport properties of the impurity atoms and show that a crossover 
from coherent to diffusive behaviour occurs with increasing interaction strength. 
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1. Introduction 

Ultracold atoms in optical lattices have attracted considerable interest during the 
last few years. Theoretical investigations showed that ultracold atoms in optical 
lattices can be used for mimicking a wide range of models encountered in condensed 
matter physics [1]. These models include the Hubbard Hamiltonian [2,3], spin- 
spin interactions [4,5], high-temperature superconductivity [6,7], effective magnetic 
fields [8-13], even with non-abelian gauge potentials [14, 15], and the fractional 
quantum Hall effect [13,16], to name but a few. Experimental efforts have led to 
an unprecedented control over the properties of optical lattice systems [17], including 
such milestones as the Mott insulator to superfluid transition [18,19], investigations of 
Bose- Fermi [20,21] as well as Bose-Bose mixtures [22], vortex pinning [23], the creation 
of repulsively bound atom pairs [24], and cold controlled collisions between atoms in 
optical lattices [25,26]. Especially for the latter experiments the absence of lattice 
phonons and thus the suppression of decoherence mechanisms was crucial. However, 
when mimicking the behaviour of electrons in crystals this lack of phonons might 
lead to an oversimplification of the underlying model and it is desirable to introduce 
phonons in a controlled manner into the optical lattice system. 

One way of achieving this goal is to immerse the optical lattice system into a 
Bose-Einstein condensate (BEG). Experiments where both atom species are trapped 
by the optical lattice are common and decoherence effects in such systems have already 
been observed [20,21]. Furthermore, these Bose-Fermi mixtures promise rich phase 
diagrams including charge and spin density wave phases [27,28], pairing of fermions 
with bosons [29] and a supersolid phase [30]. We, however, focus on the case where only 
one species is trapped by the optical lattice. This can be achieved by a suitable choice 
of the laser wavelengths and atomic species such that the BEG is not affected by the 
optical lattice [31]. The lattice atoms interact with the condensate via density-density 
interaction, which can be described by Bogoliubov phonons coupling to the impurities. 
Earlier studies on such systems have shown that this coupling can be exploited to 
cool the lattice atoms to extremely low temperatures [32,33], which are otherwise 
very difficult to achieve. In Ref. [34], the present authors have derived a model in 
which lattice atoms dressed by a coherent state of Bogoliubov phonons constitute 
polarons [35]. The model exhibits an attractive interaction potential between the 
lattice atoms [36,37] and allows a generalised master equation to be deduced which 
shows that the system exhibits a crossover from coherent to diffusive dynamics. 

Instead, in this work we concentrate on describing the lattice atoms by a quantum 
master equation (QME), which is derived in Sec. [31 We show in Sec.|3]that for the case 
of fixed impurities the system can be solved exactly and after tracing out the BEG 
degrees of freedom the QME reproduces this exact solution. Due to the coupling to 
the phonons the lattice atoms experience dephasing, which can be used to demonstrate 
the effects of sub- and superdecoherence [38] and to probe spatial properties of the 
BEG analogous to Ref. [39] . We also investigate the severe effect this decoherence has 
if the atoms in the lattice are used as a quantum register. In Sec. [S] we show that the 
attractive interaction mediated between the lattice atoms leads to the formation of 
atom clusters, which should be observable for typical experimental parameter regimes. 
This effect is reminiscent of the clustering of ad-atoms on crystal surfaces [40]. The 
influence of the phonon coupling on the transport properties of the lattice atoms is 
investigated in SeclH] We show that the crossover from coherent to diffusive transport, 
which was observed in Ref. [34], can also be described by the QME. In contrast to the 
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Figure 1. An ultracold quantum gas trapped in an optical lattice is immersed 
in a much larger BEC. The impurities interact with the BEC via the coupling 
constant k and can excite Bogoliubov phonons (right). These lead to dephasing 
effects. By exchanging phonons an off-site interaction Vi,j is mediated between 
the lattice atoms (left). 



treatment in Ref. [34], the QME in addition gives access to the off-diagonal elements 
of the density operator. We also show the limitations of the QME by applying it to a 
tilted lattice system and comparing the results to a near-exact numerical solution for 
the time evolution. 

2. Model 

The system under consideration is composed of a BEC and an ultracold gas of impurity 
atoms trapped in an optical lattice giving a setup like that shown in Fig. [TJ The 
impurities can be either bosonic or fermionic atoms. For most of the results derived 
in this work the statistics of the lattice atoms is unimportant, but for concreteness we 
will assume bosonic impurities in the following. The dynamics of the whole system is 
governed by the Hamiltonian H = Hi^ + Hb + Hi, where 



+yext(r) + |0nr)0(r) 



(1) 



describes the BEC, with (j)^ the condensate field operator, mi, the mass of the BEC 
atoms, T4xt(r) an external potential confining the BEC, and g > Q the coupling 
constant between the BEC atoms. Here, we have assumed that the optical lattice 
potential does not affect the BEC atoms, which can be achieved by choosing the laser 
wavelengths and atom species accordingly [31]. The atoms in the optical lattice, which 
are distinguishable from the atoms in the BEC, are described by the Hamiltonian iJL, 
whereas the interaction between the two sub-systems is given by 

Hi = K [ drxHr)x{r)$Hr)Hr). (2) 



Here, is the field operator of the lattice atoms, and k is the density-density coupling 
constant between the BEC and the impurities. It was shown in Ref. [34] that in 
the tight-binding limit and under the condition \K\/gno^^ <C 1 (where D is the 
dimensionality of the system, no the density of the BEC in the trap centre, and 
^ = h/ ^iribgno the healing length) the total Hamiltonian reduces to a Hubbard- 
Holstein model, given by 

Hl = -J'^a]jai + ^'^ni{ni-l) + fi'^fii, (3) 
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Hi=J2'^nuj, [Mj^X + MlXj\ fij , (4) 

V 3 

Hbec = ^1 ^^•^^IK ■ (5) 

Here, the prime at the sum indicates that zero energy modes have been excluded, 6j, 
creates a Bogohbov phonon in mode v with energy TiWj,, at creates a lattice atom in 

site j, hj = aJSj is the number operator on site j, J describes the hopping of the 
impurities, U is their on-site interaction, and the coupling between the phonons and 
the lattice atoms is given by 

Fj> = nuMj,u = K jdv (f>o{Y) [u^{y) - v^{r)] |j?j(r)|=^ . (6) 

In this equation, rjj is a Wannicr function describing an atom in lattice site j, <f>Q is 
the solution of the Gross-Pitaevskii equation and Ui, and Vi, solve the Bogoliubov- 
deGennes equations, see also Ref. [41]. 

For the case of a homogenous condensate, which we consider in the remainder of 
the paper, the expressions for the atom-phonon coupling is given by 



where q is the phonon quasi- momentum, £q = (?iq)^/2m(, the free particle 

energy, hujq = £q(eq + '^gn^) the Bogoliubov dispersion relation and /j(q) = 
dr|r7j(r)p exp(iq • r), with Q the quantisation volume. The latter integral 
can in general not be solved analytically. However, for sufficiently deep lattices [42], 
the Wannier functions can be approximated by Gaussians, yielding 

fj (q) = exp ( - 4 ^1^0 j : (8) 



with xq = ^'h/m,bUJt, where ojt is the trapping frequency of a harmonic trap 
approximating the lattice potential at a given lattice site. 

An experimental realisation of this setup is achievable with present techniques 
as follows. The creation of Rb condensates with 10^ atoms has been demonstrated 
leading to the desired BEC densities of about 10^°/m^ in three dimensions, see for 
instance [43]. By choosing a sufficiently flat trapping potential of a few Hz the BEC 
can be assumed to be homogenous to a good approximation in the centre of the trap, 
extending over a few micrometers. Furthermore, in references [20, 21] a mixture of 
Rb and K atoms has been created and trapped in a three-dimensional optical lattice, 
where several tens of thousands of lattice sites have been occupied. The filling for 
the K atoms ranged between and 1 atoms per lattice site, whereas the filling of 
the Rb atoms could exceed 5 atoms per site in the centre of the trap. Although in 
these experiments both atom species were trapped by the optical lattice, tcchniqiics 
to trap only one of the two species, for example K, have been studied extensively in 
reference [31] . Applying these techniques leaves the Rb atoms virtually unaffected and 
enables the creation of a nearly homogenous BEC. 
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3. The quantum master equation 

In order to investigate the behaviour of the impurity atoms we derive a Quantum 
Master Equation (QME) for the lattice system by tracing out the surrounding BEC. 
The details are given in |Appcndix A[ Here we note that the main condition of deriving 
the QME is that the sound velocity of the condensate, c ~ -y/gno/mt, , is larger than 
the typical hopping speed of the atoms, i.e., c » Ja/h, where a — A/2 is the distance 
between two lattice sites. With this assumption we get 

^ i ^ — A ' ^ — A sin(u.'qi) 
mdtQL{t) = Hg(t) , gL(t) -t/ / [nini^ gL[t) + gL[t)nini> - 2ni^ gL[t)ni) 

^ ^ q M' 

X {F^.iKi' + N'iiT){F^,iF*,, + F*,F^,,)) . (9) 
Here, Nq{T) — l/{exp{fLUjq/kBT) — 1) is the number of thermal phonons at 
temperature T, and ^l(^) — TrB^(t), where Tre denotes the trace over the condensate. 
Furthermore, the Hamiltonian 

H,it) ^Hl- E^E ^"sX^"'''^^ ^^'^''^'^*'' ^^^^ 
q iM 

describes the coherent evolution of the lattice atoms, most notably an off-site 
interaction which is mediated by the phonons of the BEC. The interaction term 
includes a transient behaviour described by the cosine functions, which accounts for 
suddenly turning on the interaction between lattice atoms and BEC. The sum over 
all these cosine functions vanishes in the limit t — > oo. 

The Born approximation used in deriving the QME is valid if the perturbation 
caused by the BEC is small compared to the typical energy scales given by the lattice 
system. If only the hopping term is important for the lattice atoms, this energy scale 
is determined by J. The energy scale of the perturbation is given by the so-called 
polaron energy Ep = Y^'qi^qjPqi + -^q i-Pq,0/2?i'^q ~ i^'^/g^^, which leads to the 
condition J ^ Ep. This parameter regime is complementary to the one considered in 
Ref. [34], where the opposite limit was investigated. As we will show in the following 
section, the condition J ^ Ep is not required for the case of J = 0, where an analytical 
solution of the dynamics of the whole system can be derived. After tracing out the 
BEC degrees of freedom this solution agrees with the one given by the QME. 



4. Fixed impurities and the Quantum Register 

For very deep optical lattices the hopping constant J is essentially zero and the 
impurities cannot leave their site. If there is a maximum of one atom in each 
lattice site, the setup can be used as a quantum register [44]. The lattice atoms 
represent qubits on which single qubit rotations can be implemented via external 
laser pulses. For the manipulation of the atoms single site addressability is necessary, 
which may be achieved by using infrared lattices [45] , by leaving empty sites between 
the atoms [46,47], by exploiting the properties of marker atoms [48], or by additional 
external fields [49,50]. An entangling two-qubit gate can be implemented by using 
the interaction which is mediated by the condensate [37]. For this it is necessary to 
turn the interaction on and off, which can be done by using different internal states 
of the lattice atoms, some of which couple to the BEC, and others that do not couple. 
However, due to the coupling to the BEC the qubits also experience dephasing, which 
will be investigated in this section. 



Dynamics, dephasing and clustering of impurity atoms in Bose-Einstein condensates^ 



4.I. Analytical solution of the time evolution 



We first calculate the time evolution of the atoms subject to tlie BEC coupling. For 
easier notation we focus on the case where the lattice atoms have only one internal 
state and we assume that there is a maximum of one atom per lattice site. In the 
regime were our model introduced in section [2] is valid, i.e. \K\/gnQ^^ ^ 1, and for 
J — 0, the Hamiltonian of the whole system simplifies to Hq — Hbec + Hi, where Hi 
is given in Eq. ^ and Hbec in Eq- ©■ The time evolution U{t) = cxp{~iHQt/h) of 
this Hamiltonian is solved analytically. We first note that lJ{t) — Y['^ Uq{t) with 



/7q(t) = exp 



(11) 



This operator can be decomposed using the methods described in Ref. [51], yielding 



C/q(t) = exp 



-itoqblfiqt 



exp 



. I 



X exp 



i{ujqt — aiTi{ujqt)) 



1,1' 



(12) 



Here, we have defined Fi_q — Fi^q [— isin(a;qt) — (1 — cos(a;qi))] /hujq. The last term in 



Eq. p2)) describes the off-site interaction between two atoms confined in lattice sites / 
and I' , which is mediated by a phonon with quasi- momentum q. The total interaction 
potential is derived by adding all phonon contributions together. For t oo, the 
oscillations caused by the sine function cancel each other and the off-site interaction 

3l 



potential is given b}| 



Vw = 



^'Pi,^Pi',, 



q 



(13) 



This interaction potential was already found earlier [34,37], however, the exact solution 
in addition reveals the transient behaviour after suddenly turning on the interaction, 
which is described by the sine functions. 

The coupling to the BEC leads to dephasing of the lattice atoms as gets 
apparent after tracing out the condensate degrees of freedom. There are in essence 
two different types of dephasing. One can be observed when a lattice atom is 
driven into a superposition of two different states which couple differently to the 
condensate. This situation has been discussed in reference [39]. The other type 
of dephasing is given when comparing the phase of two atoms in different lattice 
sites. We illustrate this effect of dephasing by calculating the correlation function 
{a\af3) = TT[a\apg{t)] = Tr {t)a\U{t)U^t)a|3U{t)g{0) . A combination of the 

two dephasing effects is discussed in section 

If we choose the same initial conditions as in [Appendix A[ namely g{0) = 
^b(O) and gB{0) describes a thermal state of the BEC, we get 



(al^ap) = Tri 



exp 



Hg{s)ds 



aU0)ai3{Q) exp 



Hg{s)ds] gUO) 



f For later convenience, we define this V; ;/ in such a way that the interaction term in the Hamiltonian 
is given by -J2i,i' ,l' f^lf^l' ■ 
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^'^P E' ^1^/3.1 - ^7,ql'(l - cosKi))(2iVq(T) + l)j . (14) 

Here, TrL denotes the trace over the lattice system, Hg was introduced in Eq. (|10p , and 

we made use of the identity Tre exp(Q;q6^ — aq6q)^?B ~ exp[— |Q;qp(2A'^q(T) + l)/2] 

for thermally distributed and some complex number aq [52]. The first part of 
Eq. (HH) describes the correlations between lattice sites 7 and f3 which are induced by 
the dynamics of the lattice atoms. In our case, their time evolution only contributes a 
phase term. The second part describes the dephasing which is caused by the coupling 
to the BEC atoms. Let us denote the dephasing term by T. In the thermodynamic 
limit of the Bogoliubov modes we can replace the sum by an integral and get 

r>exp(-E'(|^W) 



exp ( - / d^k ,^_^;- ^„ ' i2NUT) + 1) ) =: r^,„ , (15) 



exp(-k2a;g/2^2A 

where dq — K^rioSq exp(— q2a;§/2)/r2?ia;q, Fq.j = ■yd^exp(iqrj), and D is the 
number of spatial dimensions. For zero temperature, Z? = 3, and ^ 3> the 
integral can be approximated by neglecting the exponential function which gives 
r(T = 0) > exp(— 4K2/V27r2g2no^'^), whereas for sufficiently high finite temperature 
k^T :$> gno we find T > exp(— K^fcer/v^Trg^nQ^^). In both cases the factor T 
describing the dephasing fulfils F > <^ > for a suitable real number This is 
independent of time or of the distance between the two atoms, which is different from 
the one-dimensional case, where numerical calculations show that F decays to in 
the thermodynamic limit for f — > 00 as well as for \r-y — r[}\ 00, indicating that no 
correlations survive. 

Interestingly, the solution of the QME introduced in Sec. [3] gives the same 
correlation functions, which also holds for the time evolution of the operators Oj after 
tracing out the condensate. This is due to the fact that the interaction commutes with 
the lattice Hamiltonian for J — and that the BEC is initially in a Gaussian, namely 
a thermal, state [53]. Hence, for the case J — the QME describes the lattice atoms 
exactly after the trace over the BEC is taken. For J 7^ the interaction Hamiltonian 
and the one describing the lattice atoms do no longer commute, and the QME is only 
valid for a large hopping with J ^ Ep. 



4-2. Sub- and superdecoherence and probing the BEC 

It has been predicted in Ref. [38] that decoherence effects caused by a qubit-bath 
coupling can be enhanced (superdecoherence) or suppressed (subdecoherence) for 
certain casesQ Here we show that these effects are indeed observable in a BEC 
for realistic experimental parameters. The effect can moreover be used in order to 
probe such BEC properties as the temperature with different sensitivity, similar to 
the method introduced in Ref. [39]. 

Let us assume that the atoms have two internal states |0) and |1) which couple 
to the BEC with different coupling constants kq = 0, ki = n. This case can be 



I In our case it would be more appropriate to speak of sub- and superdephasing. We will however 
stick to the conventions introduced in Ref. [38] . 
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Figure 2. The different behaviour of the dephasing functions r_ (solid Hne), 
Fq (dashed line) and Tj^ (dotted line) versus (a) time, (b) temperature, and 
(c) distance d of the lattice atoms. Two ^^^Cs atoms were assumed in an 
optical lattice with wavelength A = 790nm and a depth of Vq = 40i?fi, where 
= (27r?i)'^/2m; A^. The one-dimensional BEC consists of ^'^Rb atoms with a 
linear number density of 5 X lO^/m. The couplings are given 

and K = 3.5 X 10~^i?RA. In (a), a temperature of T = 5nK and a distance 
of d = 5A/2 were chosen, in (b) the same distance and a time of t = 10ms, 
in (c) T = 5nK and t = 10ms. For recent experiments of RbCs-mixtures, see 
Refs. [55,56]. 



easily generalised to a non-zero kq. The different coupling strengths can for example 
be realised in a '*°K— ^^Rb mixture. The positions of the Feshbach resonances for 
the scattering with the rubidium atoms depends on the fine structure levels \F,mp) 
of the potassium atoms [54]. By choosing the external magnetic field close to one 
Feshbach resonance considerable differences between the scattering lengths and hence 
the coupling constants for different internal states can be achieved, and it is even 
possible to tune one scattering length to zero. 

Initially, the total state of the two atoms is described by the density matrix 
^2(0) — X]y7ci=o 1 9ijki{0)\'ij){kl\- After evolving for a certain time t and tracing out 
the BEC the density matrix changes to g2{t) = t^coh I]yfei=o,i ftjfe;(*)Nj)(^^l^mh' 
where Ucoh describes the coherent evolution of the atoms induced by the off-site 
interaction term, which for our case only changes the phase of the matrix elements, but 
not their absolute value. The elements 'Qijkii't) are given by 'Qijkiit) — Qijki{0) if i ~ k 
and j = /, Qi^kiit) = '^oQijkiiO) if i + j + fc + Hs an odd number, g^jkiit) = T-QijkiiO) 
for giQiio and gwoi and Qi^kiit) = 'r+gijkiiO) for gnoo and £iooii- The functions r+,o,- 
are calculated analogously to the dephasing term F in Sec. 14.11 and are given by 




j:'d. 



q 



q ^,2 \^\ 2N^iT) + l)] , (16) 



^'dq[2-2cos(q(r,-r^))]i-^^^(2iVq(T) + l)l , (17) 

q ) 

F+ = exp(-5]'dq[2 + 2cos(q(r^-r^))]i^^^(2iVq(T)-|-l)l , (18) 



where r.^ and denote the positions of the two atoms. Although the overall structure 
of these three decoherence terms looks very similar they behave quite differently. Let 
us for the moment assume that only such momenta q contribute considerably to the 
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sum for which the conditfon q(r^ — r^) <C 1 holds. This is typicahy the case for a 
temperature of T w and the atoms trapped in two neighbouring lattice sites. Then, 
r_ « 1, which means that the decoherence for an initial state |'(/') = (|10) ± |01))/\/2 
is strongly suppressed, because the fluctuations of the condensate happen on a length 
scale which is too large to resolve the distance between the two atoms. In contrast, 
the decay of the off-diagonal elements of the initial state jV') = ± |00))/\/2 is 
enhanced. These effects of sub- and superdecoherence have been predicted in Ref. [38] 
for a general type of qubit-environment coupling and are similar to the Dicke effect [57] 
well known in quantum optics. However, in our case it is not the decay rate of the 
excited state which is suppressed or enhanced. Instead, only the off-diagonal elements 
of the density matrix are affected leaving the occupation numbers unchanged. 

As depicted in Fig. [21 these effects can indeed be observed for realistic 
experimental parameters. Figure [2ja) shows the time dependence of the three 
functions. F- reaches its stationary state very quickly, whereas the two other functions 
drop exponentially to zero. The temperature dependence is shown in Fig.l^Jb). Since 
for higher temperatures phonons with shorter wavelengths become more important, 
the condition q(r^ — r^) <C 1 is no longer fulfilled and F_ decreases, however 
considerably slower than r+ . The condition is also violated when the distance between 
the two atoms is increased, as shown in Fig.[2lc). Interestingly, the function F_(_ shows 
exactly the opposite behaviour: With increasing distance, the value of the function 
increases as well, until it reaches the same value as F_. 

^.3. Implications for the quantum register 

The dephasing investigated in the previous section has severe implications on the 
performance of a quantum register. In order to store the information two internal 
states |0') and of the lattice atoms are needed, and we assume that these states 
do not couple to the BEC, in which case there is also no decoherence caused by the 
condensate coupling. However, the BEC can be used to perform a two-qubit gate 
between two atoms submerged into it as detailed in Ref. [37]. In short, if an atom is 
in state, say, it will be driven by a laser pulse to a state |1) which couples to the 
condensate, whereas the state |0') is either unaffected or driven to a state |0) which 
does not couple to the BEC. If both atoms involved are in state |1) they can exchange 
BEC phonons which causes an additional energy shift. In the previous section, this 
evolution was described by the unitary C/coh, and by appropriately chosen laser pulses 
this operator results in a controUed-phase gate. However, due to the coupling the 
atoms also get entangled with the BEC, which results in dephasing after tracing out 
the condensate. 

The influence of the dephasing on the density matrix g2 is suitably expressed by 
using Kraus operators Ej. The super-operator A giving the effect of the dephasing 
(excluding the controUed-phase gate) can be expressed as 

d 

where d is the dimensionality of the quantum system, here d = 4 for two qubits. Since 
the dephasing A commutes with the controUed-phase operation, the effect of the noisy 
gate is described by Ag(p2) = UczAi02)Ulz ~ ^0czQ-2U\z). The Kraus operators 
for A are given in |Appendix B[ where we also show how the average fidelity {F) of the 
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noisy gate can be calculated using the explicit form of these operators. We find 

(^) = ^(4 + 4ro + r_ + r+). (20) 

From the definition of the dephasing terms it is evident that the fidelity is worse for 
higher interactions k between the BEC and the lattice atoms. On the other hand, 
a lower interaction k means a lower interaction strength between the two atoms and 
hence it takes a longer time to perform the gate. In the three-dimensional case, by 
choosing an arbitrary low interaction strength k, the average gate fidelity can be 
brought arbitrarily close to I, however the time to perform the gate gets arbitrarily 
long as well. Since there also exist other decoherence mechanisms in cold atom systems, 
this will eventually decrease the performance of the setup. 

For concreteness, we consider the three-dimensional case where two ^'^■^Cs atoms 
are placed in two neighbouring lattice sites with wavelength A ~ 790nm. They 
are surrounded by a ^■^Na BEC with a number density of (5 x IO^)'^/m'^ and an 
interaction strength oi g = 1.5 x 10~'^E^_X^. Here, = (27r?i)^/2miA^ is the recoil 
energy, with mi the mass of the lattice atoms and A the wave length of the laser 
creating the lattice. The interspecies coupling is given by k = 1.1 x 10^^-ErA^. 
In the thermodynamic limit, for i ^ j and ^ 3> the mediated interaction is 
Vij = K^^exp{—^/2\i — j\a/^)/TTg^^\i — j\a, where a is the distance between two 
neighbouring lattice sites. This leads to the (minimal) gate time tg — Tni /Vi^2 = 40ms 
[37], where the gate fidelity is given by {F) = 0.99. For getting the last result, we used 
the approximations F^; > exp(— c^^K^/v^Tr^g^rio^'^), which hold in the thermodynamic 
limit, T = 0, ^ ^ xq, and where cq = 1, c_ = 2, and c+ = 4. 

For the case of two or one spatial dimensions, the interaction is in general larger 
and thus the expected gate times shorter. However, the fidelities decrease much faster 
than for three dimensions, and all of our numerical tests showed that for the same 
gate times the three-dimensional fidelities were always better than the ones achieved 
in lower dimensions. Although this restricts the applicability of this setup for quantum 
information purposes, the scheme can still be used for probing the interaction strength 
which is mediated by the condensate, simply by measuring the phase differences for 
varying interaction times. 

5. Clustering of the lattice atoms 

If the lattice is not uniformly filled with atoms but the filling ranges somewhere 
between zero and one atom per lattice site, the mediated interaction has an important 
impact on the spatial distribution of the atoms. For low enough temperatures, it will 
lead to atom clusters. To observe the clustering of the atoms we have to allow for 
a weak hopping J <C |Vi^2|, l^^lM because otherwise the atom distribution remains 
stationary. We can assume that the perturbation due to the hopping does not change 
the interaction potential Vij derived from the exact solution, which is in agreement 
with our earlier findings [34]. The hopping energy can furthermore be neglected 
compared to the interaction potential, such that the weak hopping leads only to a 
re-arrangement of the atoms, which is similar to the treatment of ad-atoms on crystal 



§ In the assumed homogeneous setup the value Vi_2 corresponds to the mediated interaction between 
two nearest neighbours and Vi,i to the mediated on-site interaction. 
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surfaces [40]. The Hamiltonian for this situation is given by 



The mediated on-site interaction Vi^i, which in our case is always attractive, 
cf. equation (fT3)) . might overcompensate the on-site interaction [/, i.e., U + Vi^i < 0. 
If this happens, all the atoms can aggregate in a single lattice site, which will inevitably 
lead to three-body losses. We therefore assume that the repulsive on-site interaction 
U is large enough such that states with more than one atom in a single lattice site can 
be neglected, i.e., U + Vi^i » J. For the parameters used below (see captions of Figs. [3] 
and [4]) the mediated on-site interaction is given by T4,i = — 0.03i?R in one dimension 
and Vi^i = — O.OS-Er in two. Since we require a deep lattice for the hopping J to 
be small, a sufficiently high U can easily be achieved. For a reasonably deep lattice 
of 15i?H and perpendicular confinement of SS-Er, the interaction strength is on the 
order of ?7 « 0.4ii^R for a one-dimensional lattice and U « 0.3£'r for two dimensions 
and thus sufficiently high to overcome the induced on-site attractive interaction. The 
repulsive interaction strength corresponds to a temperature of U /k^ — 150nK in one 
dimension {U/k^ = 130nK in two dimensions), such that neglecting states with more 
than one atom per lattice site is well-justified for the temperature regime considered, 
see below. With these assumptions, the ground state for the lattice system consists 
of a cluster where all the atoms are located in neighbouring lattice sites. We note 
that this also holds if the lattice is loaded with spin-polarised fermions, in which case 
the Pauli exclusion principle ensures that there is never more than one atom in each 
lattice site. The restriction to a maximum of one atom per lattice site and ignoring 
interactions beyond nearest neighbours makes it also possible to map the system to 
an Ising model by using the correspondence fij = (1 + Sj)/2, which gives 



where a constant term has been neglected. We will make use of this correspondence 
shortly when we compare our results to analytical ones found for the Ising model in 
two dimensions. 

In an experiment, the lattice atoms will have a finite temperature T, which leads 
to a breaking up of the ground state into smaller clusters, reflecting the increased 
average energy. It is therefore important to investigate at which temperatures the 
clustering can be observed. To this end, we simulated the lattice system using the 
well-established Metropolis algorithm [58,59], which is often used in the simulation 
of classical lattice spin and lattice gas models when the kinetic energy can be 
neglected compared to the interaction energies. For a one-dimensional lattice, the 
averaged number of cluster|j] for different temperatures is shown in Fig. [3lja). For 
temperatures above T = 7nK the number of clusters does not change very much with 
the temperature and roughly 110 clusters exist, which shows that most of the atoms 
do not have nearest neighbours. A drastic reduction of the cluster number only occurs 
for T < 7nK, such that the number of clusters finally reaches = 1 for T « 0. A 
method to achieve such low temperatures of below 7 nK has been proposed, where the 
surrounding BEC is first used for dark-state cooling [32,33]. 

The numerical data in Fig. [3] is compared to an analytical result derived in 
the thermodynamic limit, where the number of lattice sites M goes to infinity 

II Note that a single atom without nearest neighbours is also considered to constitute a cluster. 
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(22) 
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Figure 3. Normalised cluster number Nc/N versus (a) temperature of the one- 
dimensional lattice gas and (b) coupling to the BEC. The solid line shows the 
analytical estimate given in Eq. 12311 . results from numerical calculations only 
taking the nearest neighbour interaction into account are marked by 'o' and 
numerical results including the full interaction potential are marked by 'x'. The 
algorithm was run for 10* steps after which an equilibrium state is reached. Then, 
another 5 X 10^ steps were done and after every 1000 steps the number of clusters 
was calculated. Error bars show one standard deviation derived from averaging 
over those results. The lattice with wavelength A = 790nm consists of M = 800 
lattice sites filled with = 160 atoms. The surrounding BEC was assumed 
to consist of *^Rb atoms with a linear number density of n = 5 X 10^/m and 
a coupling g/E^\ = 1.1 X 10~'^. In (a), the interspecies coupling is given by 
k/E^\ = 2.5 X 10~^, whereas in (b) the temperature is fixed to T = 3nK. 
Periodic boundary conditions were assumed. 



whilst keeping the filling fraction N/M constant, and only taking nearest-neighbour 
interactions into account [60]. In this case, the normalised number of clusters N^/N 
is given by 

^ M + 4(M - N)N [exp(2|yi,2|/fcBT) - 1] /M^ - 1 
N N 2[exp(2|-t^i,2|/fcBr)-l] ' ^ ' 

where fee is Boltzmann's constant. Our calculations show that the numerical results 
taking only the nearest neighbour term of the interaction into account is in excellent 
agreement with the result for the thermodynamic limit. We furthermore observe 
that the number of clusters for the full interaction potential is for low temperatures 
considerably smaller than for the truncated one, which indicates that the interaction 
terms beyond nearest neighbours make the clusters more stable. This gets also 
apparent when the temperature is fixed and the number of clusters is calculated for 
different coupling strengths k, as shown in Fig. [3]Jb). For small k the interaction 
beyond nearest neighbours is still quite weak and does not change the number of 
clusters significantly, whereas for stronger coupling k, interaction terms beyond nearest 
neighbours are considerable, reflected in a smaller number of clusters or, equivalently, 
in a larger average cluster size. 

The situation gets more interesting if we consider a two-dimensional lattice. For 
T » 0, the atoms aggregate in an "island" and for increasing temperature parts of 
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Figure 4. (a) Distribution of the atoms in a homogenous two-dimensional 
lattice with periodic boundary conditions and a temperature of T = 1.8nK. 
A macroscopic island is clearly visible. (b) Normalised size of the largest 
cluster Ni/N versus temperature of the lattice gas for a two-dimensional system 
(marked by 'x') and normalised number of clusters Nc/N (marked by 'o'). The 
solid line shows the analytical expression given in Eq. I)24|l. the dashed line 
is a guide to the eye. Error bars show one standard deviation derived from 
averaging over several runs. The lattice with wavelength A = 790nm consists of 
M = 150 X 150 lattice sites filled with 900 "^K atoms. The surrounding BEC was 
assumed to consist of *^Rb atoms with a number density of n = 25 X lO^'^/m^ 
and a coupling g/E^,^^ = 5.1 X 10""^. The interspecies coupling is given by 
k/E^X^ = 1.87 X 10-2. 



this island break away, see Fig.[4l^a). We have calculated the size of the largest cluster 
(i.e. the number Nj of atoms constituting the island) numerically, only taking nearest 
neighbour interactions into account. This will give a lower bound for the temperature 
below which the formation of the island can be observed. The results are shown in 
Fig. Sljb). For very low temperatures, all the atoms are contained in a single cluster. 
This largest cluster gets smaller with increasing temperature and shows a pronounced 
transition at a temperature of T « 2nK, above which the system mainly consists out of 
many small clusters. It has been shown in Ref. [61], that the normalised size Nj/N of 
the island in the thermodynamic limit and only taking nearest neighbour interactions 
into accout behaves as 

iV, ^ {l + No){No-Nr) 

N 2No{l-Nr) ' ^ ' 

Here, iVo = {1 - [smh{2J,/kBT)]~^y^^, and Js = |Vi,2|/4. For this formula to be 
valid we further require a filling of N/M < 1/2, such that Nr = 1 - 2N/M > 0, 
and above the transition temperature Tj where Nj/M hits zero for the first time, the 
function has to be set to zero, since for this temperature the normalised size of the 
largest island vanishes in the thermodynamic limit. Our Monte-Carlo results are in 
excellent agreement with Eq. (I24p . deviations above the transition temperature are 
due to finite size effects. We note that a higher transition temperature T/ can be 
achieved by increasing the interaction strength Vi_j or by loading more atoms into the 

lattice. We find fee 7/ — 2Js/arsinh 1/ {/l ~ , which yields for the chosen values 
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Figure 5. (a) and (c) time evolution of tlie atom density, (b) and (d) graphical 
representation of tlie absolute values of tlie density matrix for the final time 
step t = 8.2ms shown in (a) and (c). Darker colour corresponds to a higher 
numerical value. For the examples shown we have assumed a atom in an 
optical lattice with wave length A = 790nm, and a hopping of J = 0.03-Er. The 
BEC consists of *^Rb atoms with a density of no = 5 X 10® /m and a coupling of 
g/EjiX = 1.1 X 10~^. The temperature is given by T = lOOnK. The interspecies 
coupling for case (a) and (b) was assumed to be k = 0, whereas for (c) and (d) 
it was k/E-^X = 1.94 X 10~^. For recent experiments on these mixtures, see 
Rcfs. [22,62]. 



(see caption of Fig. [?]) and iV^ — + a value of Tj = 2.34nK. 
6. Transport properties of the impurity atoms 

In this section, we use the QME to investigate the behaviour of the lattice system 
for a hopping J Ep. In general, an analytical solution of the QME is no longer 
possible, and the dynamics has to be calculated numerically. Let us consider the simple 
case of one atom initially localised in a single lattice site. With vanishing coupling 
to the condensate, the atom will spread across the lattice coherently in a wavelike 
motion, as shown in Fig. Elja). The interference fringes between the two wave fronts 
are clearly visible. The coherent nature of the evolution is also reflected in the density 
matrix, where all off-diagonal elements have their maximum possible value, compare 
Fig. [DJb) . The situation changes when the coupling to the condensate is increased, 
as shown in Fig. [5jc). The wave-fronts still exist, however the region in between no 
longer exhibits clear interference patterns, which are washed out instead. This implies 
that the coherent nature of the evolution is impaired, which is also supported by the 
density matrix, see Fig.[5ljd). The off-diagonal elements are clearly suppressed, which 
confirms the incoherent character. 

The motion of the atom in the lattice can be explained by the coexistence of 
wavelike, coherent evolution responsible for the two wave-fronts, and a diffusive, 
incoherent evolution. The diffusive motion stems from the decay of the off-site 
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Figure 6. (a) Standard deviation of the density distribution after an evolution 
time o{ t = 6.7ms (dashed line) and t = 8.2ms (solid line), (b) Average density p 
(solid line) and standard deviation p<j (dashed line). Apart from the interspecies 
coupling K, all the parameters are as in Fig. |5] 



elements in the density matrix, which destroys the memory of the atom and causes 
it to perform a random walk, leading to the washing out of the interference effects 
between the two wave fronts. To investigate the time evolution more quantitatively, 
we calculate the standard deviation of the atom density distribution at a time t, 
given by Cd = ^jPjU ~ Jo)^: where j labels the lattice sites, pj is the probability of 
finding the atom at lattice site j, and jo is the initial lattice site of the atom at t — 0. 
Figure [U^a) shows that for an increasing coupling to the condensate this standard 
deviation Ud decreases. This stems from the fact that the hopping is reduced due to 
the coupling between the lattice atoms and the phonons of the BEC, which has also 
been observed in Ref. [34] . 

Taking one standard deviation to either side of jo defines a suitable interval / 
between the two wave fronts to investigate the diffusive character of the motion. For 
this interval, we calculate the average atom density p = X^jg/ Pj /"^^a ^^id the standard 
deviation of the density distribution pd = '^j£i{Pj~py' /"^^d- For the coherent case, we 
expect that this standard deviation is on the order of the average density p due to the 
interference patterns, whereas for the incoherent case the standard deviation should 
be much lower than the average density due to the more homogenous spread of the 
atom between the two wave fronts. These expectations are confirmed by our findings 
shown in Fig.[6l[b). We observe that the average density between the two wave fronts 
p moderately increases with increasing interspecies coupling. The standard deviation 
of this density p^ stays first approximately constant with increasing k, indicating 
a regime where the evolution can still be considered to be coherent, but then for 
K > 5 X IO^^E'rA drops off considerably, caused by the increased loss of coherence. 

The QME describes the lattice atoms well as long as only dephasing has to be 
taken into account and energy exchange with the BEC can be neglected. This is 
consistent with the approximations we have used to derive the QME: It was assumed 
that Ja/h <C c, which means that the atoms move slower than the critical velocity 
of the BEC. Then, according to the Landau criterion for superfluidity, no energy 
exchange is possible with a single phonon process, and higher order phonon processes 
are not included into the QME. 

This gets especially apparent when we try to describe the decay of Bloch 
oscillations with the QME. Let us assume a Stark potential Hg = K J2j jf'-j is 
applied to the one-dimensional optical lattice. For one particle initially located in a 
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Figure 7. Bloch oscillations in an optical lattice submerged into a BEC of °' Rb 
atoms with a density of no = 5 X 10^/m, a coupling of g/E^\ = 1.1 X 10~^ and a 
temperature of T = 75nK. The atom is trapped in a lattice with wavelength 
A = 790nm and has a hopping rate of J = 0.03-Br. The strength of the Stark 
potential is given by K/h = 1.5kHz. The interspecies coupling is given by (a) 
K = and (b) k/E^\ = 1.6 X 10"^. 



single lattice site and not coupling to the BEC, we expect breathing oscillations with 
frequency lub = K/h and width 4:Ja\fiin{ujBt/2)\/ K [63]. This is shown in Fig. [Tlja). 
When the atom couples to the BEC, it can dissipate energy into the condensate 
modes and one expects that the atom can "fall down" the tilted lattice, causing a 
net current. When using the QME to describe the dynamics this is, however, not the 
case. Instead, the Bloch oscillations only get washed out, but the mean position of 
the atom remains unchanged, as shown in Fig. [7^b). This is in agreement with our 
expectations and with other attempts to describe the system with a simple Master 
Equation approach [64]. An alternative way based on a generalised master equation, 
which includes the coupling of phonons to the hopping term, has been devised to 
describe the occurrence of an atomic net current. Details will be given elsewhere [65]. 

Another way to observe the net current is given by solving the complete dynamics 
of the atom and the surrounding BEC numerically. For this, the BEC has to be 
discretised, which is equivalent to assuming that it is trapped in an optical lattice 
with appropriately chosen hopping and interaction constants [66]. The Hamiltonian 
describing the system is then given by 

Ub 



j 



(25) 



with pj (fij) an operator that creates (annihilates) a boson in lattice site j, Ub the 
on-site interaction strength and Jb the hopping matrix element of the BEC atoms, 
and Ui the interspecies interaction strength. The time evolution of this Hamiltonian 
including the BEC and a single atom in the lattice is performed numerically by using 
the time evolving block decimation (TEBD) algorithm [67-69]. The resulting density 
distribution of the lattice atom initially placed in a single lattice site and evolved for 
one Bloch period is shown in Fig. Ufa). For a weak interaction strength Ui ^ Jb the 
surrounding BEC atoms hardly affect the dynamics of the lattice atom and after one 
cycle it ends up in its original lattice site again. For increased interaction the lattice 
atom is able to dissipate more and more energy to the BEC and to drift towards lattice 
sites with a lower potential energy, which leads to a net current. This is illustrated 
in Fig. [8jb), where a competition between Bloch oscillation and net current is clearly 
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Figure 8. (a) Decay of Bloch oscillations and occurrence of a net current. Shown 
are the density distributions of one atom initially localised in lattice site 13 after 
a time evolution of t = 2iT/u!g for different interaction strengths ranging from 
Ui = {0.1; 0.5; 0.8; 1; 2} Jf, (top to bottom line in lattice site 13). The dynamics 
of the lattice atom and the BEC consisting of 25 atoms was calculated in an 
optical lattice with 25 lattice sites using the TEBD algorithm. We used Ja = Jbt 
Ut = 0.5J(,, K = 1.25Jf,, T = 0. (b) Time evolution of the density distribution 
for one lattice atom with Ui = Ji,. The coexistence of Bloch oscillations and a 
current is visible. The vertical line indicates the time where the densities of (a) 
are taken. 



visible. The current can be measured either directly or after a time of flight expansion 
of the lattice atoms as detailed in Ref . [70] . 

7. Conclusion 

In the present paper we have investigated the behaviour of an optical lattice system 
immersed into a BEC. To this end, we have derived a quantum master equation which 
describes the time evolution of the lattice atoms. For the case of fixed impurities, 
the lattice system represents a quantum register and the internal states of the atoms 
represent the states of the qubit. We have derived an exact solution for this system 
which is reproduced by the QME. We found that an interaction between the lattice 
atoms is mediated by the condensate and that the coupling to the phonon modes of the 
BEC causes a dephasing of the lattice atoms, i.e., of the qubits they represent. This 
dephasing can be used to probe the properties of the BEC with different sensitivity. 

The interaction between the lattice atoms leads to a clustering process at low 
enough temperatures. In contrast to on-site clusters, which are prone to three-body 
losses, the clustering in our case is caused by the off-site interaction terms. We 
have simulated the clustering process using the Metropolis algorithm and compared 
our results with analytical ones originally found for spin systems or ad-atoms on a 
crystal surface. Our findings indicate that for realistic experimental parameters the 
clustering process occurs at temperatures of about 5nK and is observable with near- 
future experimental techniques. 

We also investigated the transport properties of the lattice atoms subject to the 
BEC coupling. We found that the lattice atoms undergo a crossover from coherent 
to incoherent evolution when the interaction to the condensate is increased. This 
process is indicated by a washing out of the typical interference fringes in the density 
distribution and by vanishing off-diagonal elements in the one-particle density matrix. 
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However, due to the approximations necessary to derive the QME, it does not include 
dissipation of energy to the BEC. This process is a vital effect to describe the decay 
of Bloch oscillations in a tilted lattice. More work is necessary to include this effect 
into a master equation description. 
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Appendix A. Derivation of the QME 

The details of deriving the QME are given in this Appendix. Our starting point is the 
Liouville-von Neumann equation for the total density operator Q{t), 

mdtm = [H,Q{t)]- (A.i) 

For a general operator O in the Schrodinger picture we define the corresponding 
operator 0{t) in the interaction picture by 

d{t) = S"B+HL)t/h Q ^-i(Hu+H^)t/n ^^ 2) 

With this, the Liouville-von Neumann equation changes to 

ihdtm^[Hi{t)rQ{t)]. (A.3) 

If there was no interaction between the lattice atoms and the BEC before a time 
to — 0, it is reasonable to assume that the density operator at time t^ is described by 
^(0) = ^l(O)® ^b(O), i.e., at time t^ there are no correlations present between the BEC 
and the lattice atoms. Here, ^l(O) and ^b(O) are the initial density operators of the 
lattice system and the BEC, respectively. We furthermore assume that at to the BEC 
is in a thermal state with temperature T . Integrating the Liouville-von Neumann 
equation, substituting it into itself and using the definition of Hj(t), we find after 
applying the Born approximation and taking the trace over the BEC variables [52] 

dtgL{t) = --2 / d^' Y.^Mt)ni'it')QL{t') - m,it')gLit')hiit)]C^^u'it - t') 

+ [gL{t')ni,it')Mt) - m{t)gL{t')m,it')]C*^^,^,{t ~ t')} , (A4) 

with 

Cq,M' {t ~ t') = F^jF:^ i, (1 + Aq(r)) e-''^-^*-*') + F:^jF^j,N^{T) e'--(*-*') . (A.5) 



Under the assumption that the sound velocity of the condensate, c ^ \J grif^lrn^, , is 
larger than the typical hopping speed of the atoms, c ^ Ja/Ti, the dynamics of the 
BEC is much faster than the typical dynamics of the lattice atoms and we perform 
the Markov-approximation by replacing gL{t') — > gh^t) and ni{t') — > fiiit). After 
calculating the time integral and transforming back into the Schrodinger picture we 
finally arrive at Eq. 
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Appendix B. Kraus operators for the two-qubit dephasing 

Here, we give explicit formulas for the Kraus operators describing the operator A given 
by A{g2) = ulzMMUcZ. We find 



El 


= + 


2ro + r+)/4 14, 


(B.l) 


E2 


= V(i- 


2ro + r+)/4(T^®CT^, 


(B.2) 


E3 


= V(i- 


r_)/4a,® I2, 


(B.3) 


Ei 


= V(i- 


r^)/4i2®(T, , 


(B.4) 




= v/(r- 


-r+)/4diag(-l, 1,1,1), 


(B.5) 


Ee 


= V(r- 


-r+)/4diag(l, 1,1,-1), 


(B.6) 



and all other Kraus operators are zero. Here, 1„ is the n x n unit matrix, cr^ is the 
Pauli z matrix, and diag is a diagonal matrix with the entries on the diagonal given 
as the argument. 

With the knowledge of the Kraus operators the average fidelity (F) of the gate is 
calculated as follows. We first note that the fidelity F of the noisy operation Ag with 
respect to the perfect (i.e., noise- free) operation Ucz on a special input state is 
defined by 

By taking the average over all possible input states the overall performance of the 
gate operation can be estimated. It has been shown [71,72] that the average fidelity 
is given by 

<^> = i,-. ^<*> = 5(dTT) (I I "<^''l^ * ") • 

where the integral is taken over the unit sphere S'^'^^'^ embedded in 2d — 1-dimensional 
real space, which is isomorphic to the d-dimensional complex space after eliminating 
a global phase, and dip is the normalised measure over the sphere, also know as 
Haar measure. With this formula and the explicit form of the Kraus operators the 
calculation of the average fidelity Eq. (PO)) is straight forward. 

We note that the Kraus operators for a coupling to independent reservoirs are 
given by replacing r_ = r+ = (Fq)^. This gives an average fidelity of (i^ind) = 
(2 + 2ro + (ro)^)/5. Hence the coupling to the BEC reservoir is worse than the one 
to independent reservoirs if r+ + r_ < 2(ro)^. 
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